G3 Genes|Genomes|Genetics

The evolutionary mechanisms leading to nonfunctional trait loss are not well understood. Some Arabidopsis thaliana populations display incomplete loss of short stamens, floral organs that do not contribute significantly to seed set in this species. In nature, short stamen number is correlated with ovule number and flowering time, suggesting that pleiotropy with these traits could drive stamen number evolution. To investigate the role of pleiotropy in short stamen number evolution, we performed QTL mapping of stamen number and ovule number and examined previously published QTL analyses of flowering time for two sets of RILs, Belm-12 x Roda-47 and Tsu-1 x Kas-1. Flowering time, ovule number, and short stamen loss were correlated in Belm-12 x Roda-47 RILs, but not Tsu-1 x Kas-1 RILs. For both stamen and ovule number, some QTLs were unique to each RIL set and others were shared. Notably, for Belm-12 x Roda-47 RILs, we mapped QTLs affecting ovule and short stamen number to a region of chromosome 5 that also affects flowering time. For Tsu-1 x Kas-1 RILs, we identified the same chromosome 5 region as a short stamen number QTL, but ovule number and flowering time QTLs map to other genomic regions. While some genetic architectures of short stamen number and traits that are correlated in nature suggest linked or pleiotropic loci on chromosome 5, other architectures allow for recombination between QTLs that affect stamen number, ovule number, and flowering time. The role of pleiotropy in short stamen number evolution is therefore context dependent.

UNC-6/Netrin is a conserved regulator of dorsal-ventral axon and cell migrations. UNC-6 is composed of a Laminin N-terminal domain (LN), three epidermal growth factor repeats (EGF), and a Netrin C terminal domain (NC). Here, we identified missense mutations in distinct UNC-6 domains and assessed their roles in dorsal VD/DD motor axon guidance and ventral AVM axon guidance. A missense mutation in a conserved residue of the LN domain (G289D) resulted in dorsal and ventral axon guidance defects similar to unc-6 null. A distinct missense mutation in the LN domain (S120F) was hypomorphic and strongly perturbed ventral AVM axon guidance with minimal effects on dorsal VD/DD axon guidance, showing that S120F is predominantly required for ventral guidance. Missense mutations altering conserved cysteine residues involved in di-sulfide bonding in the EGF domains were analyzed. EGF1(C321G) caused both ventral and dorsal axon guidance defects albeit weaker than unc-6 null, indicating that EGF1 is required for both. EGF2(C347Y) strongly affected dorsal VD/DD axon guidance similar to unc-6 null, with weaker perturbation of ventral AVM axon guidance. Previous results revealed that EGF3(C410Y) specifically disrupted dorsal axon guidance, a result that we confirmed. Our studies using missense mutations in the endogenous unc-6 locus complement previous structure-function studies using transgenic expression, and identify domains specifically required for ventral AVM guidance (S120Y in the LN domain) and dorsal VD/DD axon guidance (C410Y in EGF3). The crystal structure of UNC-6 indicates conserved N-linked glycosylation at N114 and N128. Mutation of these sites in UNC-6 had no effect on dorsal ventral axon guidance, showing that they do not play a major role. However, the N114 and N128 mutations interacted genetically with unc-40 and unc-5 mutations, indicating that these glycosylation sites indeed have a role in UNC-6 signaling. Our results will inform studies on how these distinct UNC-6 domains interact with guidance receptors (e.g. UNC-40/DCC and UNC-5) and other extracellular molecules to mediate dorsal-ventral axon guidance.

The Tyr::CreERT2 transgenic mouse lines are essential tools for conditional gene manipulation in melanocytes and are widely used in melanoma research. Two independent lines are in common use: the Bosenberg line and the Larue line. Precise knowledge of transgene integration sites is critical for designing complex genetic crosses, yet the integration sites for these lines have only recently been characterised by whole genome sequencing. Here we report that a Tyr::CreERT2 mouse stock routinely used in BrafCA;Ptenflox melanoma models carries a previously undescribed integration on Chromosome 1, distinct from the previously reported Chromosome 2 integration. Using long-read nanopore sequencing, we mapped the integration to an intergenic locus between Alppl2 and Alpi, revealing a 2,430 bp genomic deletion at the insertion site. We developed position-specific junction PCR and qPCR assays to genotype this allele and confirmed that offspring are born at Mendelian ratios. Given that this stock is associated with elevated spontaneous melanoma penetrance, accurate characterisation of this allele has direct implications for the many laboratories employing this widely distributed melanoma model. SIGNIFICANCEThis study identifies a previously undescribed Tyr::CreERT2 transgene integration on Chromosome 1 in JAX strain 013590, a line widely used across the melanoma research community. Accurate characterisation of this allele has direct practical implications for the many laboratories that rely on this model for studies of melanoma initiation and progression.

Gene model for the ortholog of raptor in the D. grimshawi May 2011 (Agencourt dgri_caf1/DgriCAF1) Genome Assembly (GenBank Accession: GCA_000005155.1) of Drosophila grimshawi. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Genome-wide association studies (GWAS) use statistical models to correlate single nucleotide polymorphisms (SNPs) to a phenotype of interest. This scan of the entire genome identifies regions of association with a phenotype, but due to linkage disequilibrium (LD), GWAS on their own cannot identify single genes responsible for phenotypic variation. Rather, fine-mapping of GWAS regions is required, necessitating the use of additional tools and software. With the introduction of more pangenomic resources in a number of crops (Guo et al. 2025; Hufford et al. 2021), the fidelity of these fine-mapping efforts is growing, presenting the opportunity to leverage new information about allelic variation towards gene discovery (Shi et al. 2023; Della Coletta et al. 2021). Panvar is a tool developed to integrate existing software and resources to perform GWAS and fine-mapping in one seamless step. For each identified GWAS peak, panvaR outputs information about LD and SNP effect prediction for each SNP and by layering locations of nearby genes, creates a refined list of possible candidate genes. We have implemented Panvar as an R package, "panvaR", which runs the analysis functions, creates interactive and static visualizations, and outputs results tables. This tool seeks to bridge the gap between GWAS and gene speeding up an important step of quantitative genetic studies.

Gene model for the ortholog of raptor in the May 2011 (Agencourt Dere_CAF1/DereCAF1) Genome Assembly (GenBank Accession: GCA_000005135.1) of Drosophila erecta. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of tango (tgo) in the May 2011 (Broad dper_caf1/DperCAF1) Genome Assembly (GenBank Accession: GCA_000005195.1) of Drosophila persimilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of Density regulated protein (DENR) in the Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of D. eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Gene model for the ortholog of tango (tgo) in the Sep. 2015 (UC Berkeley ASM127793v1/DbusGB1) Genome Assembly (GenBank Accession: GCA_001277935.1) of Drosophila busckii. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Lima bean (Phaseolus lunatus) is a broadly adapted, economically important leguminous crop and a susceptible host of root-knot nematodes (Meloidogyne spp.; RKN), which are a devastating plant pathogen in agricultural systems worldwide. To date, there have been few studies to elucidate the genetic determinants of RKN resistance in lima beans. Understanding the genetic mechanisms underlying resistance is essential for improving resistance traits and incorporating them into lima bean breeding programs. To assist in marker-assisted selection, we aimed to identify and map quantitative trait loci (QTLs) conferring RKN resistance-related traits. Three recombinant inbred line (RIL) populations were used in this study. Three populations were derived by crossing two RKN-resistant parents with the same RKN-susceptible parent and with each other. All populations were genotyped using genome-wide single-nucleotide polymorphism (SNP) markers. Each population was screened for root galling (RG) and RKN egg reproduction (ER) in response to M. incognita and M. javanica in greenhouse experiments. Three major QTLs were detected and mapped on chromosome Pl04 (QRk-pl04.1), Pl05 (QRk-pl05.1) and Pl10 (QRk-pl10.1) across populations. Among them, QRk-pl05.1 and QRk-pl10.1 affected levels of RG and ER of both RKN species, while QRk-pl04.1 suppressed root galling and reproduction responses of M. incognita but not of M. javanica. These chromosomal regions defined by flanking markers will help guide marker-assisted breeding and gene discovery for broad-based RKN resistance in lima beans.

Advances and accessibility to next-generation technology provide opportunities to sequence non-model organisms. Despite this increase in whole-genome sequencing data, annotation and the production of reference genomes remain limited. Reference genomes are a critical tool for a variety of studies in evolutionary biology, functional genomics, and conservation genetics. Tiger beetles (Cicindelidae) are a diverse and globally distributed family of beetles that serve as bioindicator taxa and flagship species for insect conservation. Here, we report highly complete, contiguous, and annotated genome assemblies representing draft reference genomes for three species of tiger beetle spanning the phylogeny. These draft reference genomes are for Audouins night-stalking tiger beetle, Omus audouini; the montane giant tiger beetle, Amblycheila baroni; and the western red-bellied tiger beetle, Cicindelidia sedecimpunctata. Article SummaryTiger beetles are a charismatic group with [~]3000 species distributed globally. Despite their popularity among insect enthusiasts and their role as bioindicators of ecosystem health, the group currently lacks a reference genome. This article outlines genome assembly and annotation for three tiger beetle species that span evolutionary relationships within the lineage. Quality control analyses show that the assemblies are reference quality and demonstrate high contiguity, completeness, and accuracy. The resulting draft genome annotations will be a valuable resource for scientific endeavors and allow for continued research on tiger beetles and their allies.

Haplotype-based breeding strategies have emerged as promising approaches to maximize long-term genetic gain by identifying complementary parental combinations while maintaining genetic diversity. However, these methods typically require phased genotypes and more intensive workflow pipelines and skillsets. We developed a novel local genomic estimated breeding value (localGEBV) fitness function with similar intent to the optimal haplotype stacking (OHS) framework fitness function and implemented both in the novel R package, HapSelect. Our aim was to evaluate whether phased haplotypes provide additional benefit over the more easily available dosage-based unphased genotypes in highly inbred crops. A subset of bread wheat nested association mapping (NAM) population comprising 444 lines genotyped with 6,054 DArT-Seq markers was analysed. Marker effects were estimated using rrBLUP, localGEBV and haplotype effects were calculated across linkage disequilibrium-defined haploblocks, and genetic algorithms (GA) were used to identify optimal sets of 30 founders using either a localGEBV derived fitness function with unphased, dosage inputs or the OHS fitness function with phased inputs. Selected parental sets were compared with conventional truncation selection (TS) through 150 generations of forward simulation. The OHS fitness function achieved a marginally greater optimized ultimate GEBV than the localGEBV fitness function during GA optimization, with only 18 of the 30 selected founders overlapped between the two methods. Despite these differences, forward simulations demonstrated nearly identical long-term genetic gain for localGEBV and OHS-selected founders, with both approaches outperforming conventional truncation selection by maintaining greater genetic diversity and delaying the genetic plateau. The minimal difference between localGEBV and OHS is likely attributable to the high homozygosity of the population, where localGEBV and haplotype effects are nearly confounded. These results demonstrate that dosage-based localGEBV provides a practical alternative to phased haplotype approaches for parent selection in inbred crops, substantially simplifying genomic workflows while maintaining long-term breeding performance. Future work should evaluate these methods in more diverse inbred populations and outbred species, where great haplotypic diversity may increase the advantage of true haplotype-based optimizations.

The Candida Genome Database (CGD; www.candidagenome.org) is both a model organism database and a fungal pathogen database. As a model organism database, CGD stores data for Candida albicans, which serves as a model species both for other Candida spp. and for non-Candida fungi that form biofilms and undergo routine morphogenic switching. As a fungal pathogen database, CGD now hosts locus pages for six species of the best-studied pathogenic fungi in the Candida group. Pathogenic Candida species have become increasingly drug resistant and there is thus a pressing need for research into basic Candida biology, epidemiology, phylogeny, and potential new antifungals, as well as a single location where all of the available data are collected, curated, and made easily searchable. CGD curates the gene-based Candida experimental literature in real time, extracting, organizing and standardizing gene annotations. CGD also links clinical data on disease to relevant Literature Topics to improve searchability for clinical researchers. Because CGD curates the literature for multiple species and most research focuses on aspects related to pathogenicity, we focus our curation efforts on assigning Literature Topic tags, collecting detailed mutant phenotype data, and assigning controlled Gene Ontology terms with accompanying evidence codes. Our Summary pages for each locus include the primary name and all aliases for that locus, a description of the gene and/or gene product, detailed ortholog information with links, a synteny view, a JBrowse window with a visual view of the gene on its chromosome, links to Phenotype, Gene Ontology, Interactions, and Expression pages, as well as sequence information, references cited on the summary page itself, and any locus notes. The database also serves as a community hub, where we link to various types of reference material of relevance to Candida researchers, including colleague information, news, and notice of upcoming meetings. We routinely survey the community to learn how the field is evolving and how needs may have changed. Here we describe CGDs new modern web interface and multiple new tools that have been added in the last 6 months, allowing, among other things, users to better understand the available expression data for a locus and seamlessly switch between species for a given locus.

Anthropogenic copper (Cu) contamination from agriculture, mining, and industrial runoff creates environmental gradients affecting physiology and behavior in wild populations. While Cu toxicity in Drosophila melanogaster is well characterized, it remains unclear whether Cu resistance is one integrated trait or several independently evolving components. Using a subset of recombinant inbred lines (RILs) from the Drosophila Synthetic Population Resource (DSPR), we measured three components of Cu response: feeding avoidance, oviposition avoidance, and physiological tolerance (median lethal time, LT50) under sustained Cu exposure. All three traits showed substantial phenotypic variation among RILs. Feeding and oviposition avoidance were both highly heritable (H 2 ~ 0.88), and RIL identity accounted for 49.5% of the variance in LT50. However, the three traits showed no significant correlation across RILs, indicating distinct genetic architecture. We identified a single male specific quantitative trait locus (QTL) on chromosome 2R that explained 17.7% of the variation in feeding preference; the interval included candidate detoxification genes Jheh1, Jheh2, Jheh3 and sano, the latter of which is associated with olfactory behavior. No significant QTL were detected for oviposition preference, suggesting a highly polygenic structure that may difficult to detect with our limited panel size. Together, these results indicate that Cu resistance in D. melanogaster is genetically modular. Behavioral avoidance during feeding, oviposition, and physiological tolerance are heritable but architecturally distinct components, each with potential to respond to selection independently.

Testing for distinctness, uniformity, and stability (DUS) is a requirement for plant variety registration and based on phenotypic traits, which is time-consuming and sensitive to environmental variation. Advances in genomics allow to complement DUS testing with molecular markers, for which two models in DUS testing were proposed by the Union for the Protection of New Varieties of Plants (UPOV). A use cases was described for maize, but an implementation has been hindered by a lack of suitable markers and validated analytical frameworks. We address these challenges by integrating historical DUS characteristics scores from 352 European hybrid maize varieties with high-density genome-wide single nucleotide polymorphism (SNP) data. Using genome-wide association studies (GWAS), we identified 18 genomic regions and candidate genes associated with 12 DUS characteristics, enabling the development of diagnostic markers consistent with the UPOV model "Characteristic-Specific Molecular Markers". Since most DUS traits are polygenic, we combined GWAS-informed marker selection with XG-Boost-based machine learning to predict notes of DUS characteristics. This approach achieved strong predictive performance across multiple traits (mean accuracy 0.67), demonstrating its potential for managing reference collections under UPOV model "Combining phenotypic and molecular distances in the management of variety collections". Both approaches were validated for two characteristics using independent public USDA-NPGS maize datasets (>1,700 accessions) highlighting the value of public data for method validation. We also identify key limitations of historical DUS data, including imbalanced and sparse trait representation, and discuss mitigation strategies. Despite these constraints, our results demonstrate that molecular markers may improve maize DUS testing, enabling faster, more accurate variety registration and supporting accelerated crop improvement. Key messageHistorical DUS datasets can be used to identify marker-trait associations of DUS characteristics using genome-wide association study (GWAS) and to develop a genomic prediction framework for an accurate prediction of DUS character notes from marker data.

Degenerate sex chromosomes (e.g., the Y or W) remain among the most difficult regions of eukaryotic genomes to assemble because they are highly repetitive and structurally complex. While many recent lepidopteran genome assemblies contain W chromosome scaffolds, the accuracy and consistency of these assemblies remain uncertain, due to lack of replication within species. However, the silkworm moth Bombyx mori is an exception, with numerous independent W chromosome assemblies currently available. We compared six independent long-read W chromosome assemblies, which proved to be highly inconsistent in structure, even among nominally identical genotypes. In contrast, autosomes and the Z chromosome were highly concordant among these assemblies, indicating that current assemblies remain unreliable for resolving W chromosome structure. Additionally, we analyzed repetitive DNA content across the genome. First, we combined assembly- and read-based repeat-discovery methods to generate a comprehensive and curated Bombyx repeat library, which we make publicly available. Assessing repeat content and diversity, we find that the W chromosome is comprised almost entirely of repetitive DNA but that the richness and divergence of W repeats are substantially reduced compared to the remainder of the genome. This reduced diversity, initially inferred from assemblies, is confirmed by direct analysis of PacBio HiFi sequencing reads partitioned by chromosome. We also demonstrate that the B. mori p50ma genome assembly (the current NCBI RefSeq assembly) carries a W chromosome and mitochondrial genome introgressed from B. mandarina. This discovery provided an opportunity to investigate patterns of divergence between closely related W haplotypes, revealing substantially more rapid turnover of repeat content on the W than elsewhere in the genome. Together, our results show that current W chromosome assemblies, although structurally flawed, nevertheless capture robust biological patterns of repeat diversity and support the hypothesis that rapid repeat turnover, rather than frequent chromosome replacement, may underlie the apparent lack of W chromosome homology across Lepidoptera.

During meiosis, accurate chromosome segregation requires significant condensation and compaction. These processes are mediated by condensins, cohesins, and histone tail modifications. We identified that MET-2, a histone methyltransferase that catalyzes the dimethylation of histone H3 lysine 9 (H3K9me2), differentially impacts chromosome size in the male vs. female C. elegans germline. In met-2 null worms, autosomes during spermatogenesis are significantly larger than wild type, while chromosome size during oogenesis is unaffected. X-univalent size in males is also unaffected by loss of MET-2, indicating MET-2 differentially regulates autosomal and X-chromosome compaction in male spermatogenesis. Autosome size is not changed when males harbor a catalytically deficient MET-2 (met-2CD) or have mutations preventing germline histone H3K9 methylation (H3K9R). In addition, met-2 males, in contrast to met-2CD or H3K9R males, have more active RNA pol II in later stages of meiosis. These data suggest MET-2 plays a noncatalytic role in mediating chromosome structure and transcription. In met-2 male germ lines, genes on the X chromosome, which is typically enriched in H3K9me2, are significantly more likely to be upregulated than genes on autosomes, even though X-univalent size is unchanged. These results suggest that MET-2 plays a sex-specific role that is not limited to its enzymatic activity.

Gene model for the ortholog of raptor in the D. eugracilis Apr. 2013 (BCM-HGSC/Deug_2.0) (DeugGB2) Genome Assembly (GenBank Accession: GCA_000236325.2) of Drosophila eugracilis. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Receptor tyrosine kinases (RTKs) are cell surface proteins that govern many critical cell fate decisions and their dysregulation is a major cause of diseases such as cancer. Much of what we know about how these proteins work in cells and tissues comes from model organisms such as the fruit fly Drosophila. Here, we identify and characterise a previously unstudied Drosophila receptor tyrosine kinase encoded by CG3277, which we name beanbag (beba). Ectopic beba expression activated Akt and ERK phosphorylation and produced gain-of-function phenotypes resembling those caused by other Drosophila RTKs. Using a MiMIC-derived T2A-GAL4 allele, we show that Drosophila beba is expressed in digestive, nervous and reproductive systems, in locations suggestive of potential roles in endoreplication and/or stem cell niche support. Animals transheterozygous for beba loss-of-function alleles were viable, developed at a normal rate, and showed no detectable change in enterocyte DNA content under standard conditions. However, beba loss-of-function females had fewer ovarioles, consistent with a role in the ovarian terminal filament, and males had increased testis hub cell number and hub volume, suggesting beba may regulate somatic niche architecture in the Drosophila gonad. Phylogenetic analysis places Beba within a Ret/Tor-related RTK radiation and supports the existence of a distinct Beba family in insects. Together, our data define Beba as a lineage-restricted Drosophila RTK with specialised roles in reproductive niche organisation.

Acrylamide is a probable dietary carcinogen formed through high-temperature cooking processes from one of its precursors, free asparagine (Asn). To reduce acrylamide formation potential in maize-based food products, the genetic basis for the accumulation of free Asn in maize kernels must be elucidated. In this study we integrated complementary quantitative genetic approaches including quantitative trait locus (QTL) mapping of a biparental population and genome-wide association studies (GWAS) in a diverse inbred panel. We ultimately identified one major QTL for free Asn which was not reflected in the GWAS results. Subsequent sequence analysis of the QTL mapping population parents revealed a 921bp deletion in As-paragine Synthetase 3 (ASN3) in the low Asn parent. Validation is underway with rtPCR, genotyping the association panel for the deletion, and evaluation of a CRISPR-Cas9 knockout of ASN3. This work provides a foundation for breeding and genome-editing strategies to reduce acrylamide-forming potential in maize-based foods.